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The average absorption properties of broad absorption line 
quasars at 800 < Arest < 3000 A, and the underlying physical 
parameters 
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ABSTRACT 

Broad absorption line quasars (BALQs) present a large diversity in their broad ab- 
sorption line (BAL) profiles. To investigate what physical parameters underlie this 
diversity, we construct a sample of BALQs which covers Arcst ~ 800 — 3000 A, based 
on the Sloan Digital Sky Survey DR7 quasar catalogue. The average BAL properties 
are evaluated by taking the ratios of average BALQ spectra to the average spectra 
of matched samples of non-BALQs, where the matching is based on various emission 
properties. We find the following properties, (i) There is no detectable Lyman edge 
associated with the BAL absorbing gas (r < 0.1). (ii) The known increase of average 
absorption depth with the ionization potential extends to the higher ionization Nv 
and O VI BALs. We also find that the C iv BAL profile is controlled by two parameters, 
(i) The Hell emission EW, which controls the typical velocity of the Civ BAL, but 
does not affect the absorption depth, (ii) The spectral slope in the 1700-3000 A range 
(auvi), which controls the C IV peak absorption depth, but does not affect the typical 
velocity. The He ii EW and auvi also control the observed fraction of quasars that are 
BALQs. We suggest that a lower Hen EW may indicate a weaker ionizing continuum, 
which allows the outflow to reach higher velocities before being overionized, possibly 
without a need to invoke a shielding gas. A redder continuum may indicate a more 
inclined system, and a larger covering factor and larger column of the outflow along 
the line of sight. 

Key words: galaxies: active - quasars: absorption lines - quasars: general. 



INTRODUCTION 



Broad absorption line quasars (BALQs) are a subtype 
of active galactic nuclei (AGN), defined by the pres- 
ence of broad and blue-shifted ab sorption features (e.g., 
IWevmann. Carswell fc SmithlligSlh . The intrinsic fraction 
of BALQs from the total quasa r popu lation is csti- 



mated to be ~ 15 — 20 per cent (IHewett fc Folta 



2002; 



200£. 



Reichard et al.ll2003l : JKnigge et al.|[2008l : iGibson et al.l 

Allen et al.l I2OIII claim it can be as high as 

per cent). While there are several differences i 

sion properties between BALQs and non-BALQs, BALQs 

appear to be drawn from the non-BALQ population 



- 40 
emis- 



high-ionization lines (these objects are termed 'HiBALQs'). 
For a smaller fraction of BALQs, absorption is also ob- 
served in low- ionization lines, e.g. Mgll ('LoBALQs'). The 
predominant unifying model states that the difference be- 
tween the AGN subtypes is our viewing angle towards the 
quasar central regions (e.g., [Elvig [2000|). However, an al- 
ternative scenario suggests that LoBALQs are an evolu- 
tionary stage of AGN, in which the nucleus expels a sur- 
rounding dusty 'cocoon' (IVoit. W evmann fc Korist al Il993l: 
Urrutia et alll2009l: iFarrah et al.ll2010 . 2012 : Glikman et all 



20121 : cf iLazarova et al.ll2012l ') 



JWevmann et al.l Il99ll : iHamann. Korista fc Morris! 1 19931 : 
iReichard et al.ll2003h . Broad absorption lines (BALs), and 
the Civ BAL in particular, span a large range in depth, 
width and overall velocity shift («shift) between different 
objects. For most BALQs, absorption is observed only in 

* E-mail: alexei@physics.technion.ac.il 



What are the average BALQ spectral properties short- 
ward of LyQ, and near the Lyman limit? Most of 
BAL Q studies analyse spectra only down t o Arest of Siiv 
(e.g.. [Gibson et al.ri2009( : lAUen et all I2OIII ) or Nv (e.g., 
IWevmann et al.lll991l ). In t his study, we uti lize the Sloan 
Digital Sky Survey (SDSS; York et al. '2OOO') DR7 quasar 
catalogue (^Schneider et al. 2010 ; Shcn et al. 2011), and in- 
vestigate high-2 (> 3) object spectra, which cover A,cst from 
C IV down to ~ 800 A. The spectra of z > 3 quasars 
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are heavily absorbed by the intervening Lyman forest at 
Arcst < 1216 A. We overcome this foreground absorption and 
derive the BALQ intrinsic absorption by taking the ratio of 
the average spectrum of BALQs and the average spectrum 
of a matched control sample of non-BALQs. The ratio spec- 
trum allows us to place a limit on the average Lyman edge 
depth associated with the BAL systems, and thus a limit 
on the covering factor (CF) of low ionization BAL system. 
We also extend earlier studies on the average absorption 
strength of C ll, Si iv and C iv to the higher ionization N v 
and Ovi BALs. 

The ratio spectrum is meaningful only if the aver- 
age emission properties of BALQs and non-BALQs are in- 
deed the same. Various studies find systematic differences 
between the intrinsic emission properties of BALQs and 
non-BALQs. There are reports that BAL Qs are located on 
the hi gh-L and high-L/ I/Edd end of th e iBoroson fc GreerJ 
l|l992ll eigenvector 1 (Boroson| bOOd ). BALQs are ob- 
served to have larger blueshift s of Civ emission than non- 
BALQs (JRichards et al.l |2002| : iReichard et"la] 120031), and 
LoBALQs to have the highest blueshifts I Richards et al.l 
120021). BALQs are fo u nd to be redder tha n non-BALQs 
(IReichar d et al.1 l2003l: iMaddox erall l2008l : iGibson et al] 
[2009 ; AUeri et aLlbOHl). aiid LoBALQs to be redder than 
HiBALQs llWevmann et al.lll99ll: ISpravberrv fc Foltj[l993: 



Boroson fc Meverslll992l : lReichard et al.ll2003l : iGibson et all 
20091 ). The reddening is interpreted by some authors as a 



result of BALQs being vie wed preferentially closer to edge 
on (e.g., I Ogle et al.l Il999f ). or a s an imprint of a dusty 
'cocoon' (see references above). iTrump et al.l (|2006l ) find 
the emission lines to be broader for BALQs than for non- 
BALQs. BALQs are reported to have lower Civ equivalent 
width (EW) than non-BALQs, X-ray weaker BALQs are 
observed to have stronger BALs with larger terminal ve- 
locities, and t he measured velocities are l a rger for higher- 
Luv BALQs ("Brandt. Laor fc Wills! |200G| : iLaor fc Brandt 
l2002l : lGibson ct al. 200 j). LoBALQs sho w the strongest ari d 
broadest high-ionization absorption lines iJAUen et al.ll201ll ) . 

Thus, to make a more accurate ratio spectrum one needs 
to use samples of BALQs and non-BALQs with a similar dis- 
tribution of intrinsic emission properties. One can then take 
another step, and explore whether the derived absorption 
properties, in particular the average C iv absorption profile, 
depend on the intrinsic emission properties, such as L/L^dd- 
For this purpose, we expand our study to z ~ 1.5 quasars, 
where one can observe the 1400-3000A range, and derive 
emission parameters, such as L/L^dd, based on the Mgll 
A2798 broad emission line. Despite numerous studies, the 
intrinsic properties which underlie the diversity of BALs re- 
main elusive. This study allows us to address what causes 
the large diversity of the observed BAL properties 

We explore in this study the dependence of the C iv 
BAL properties on the Hell EW and on the UV slope 
(auv). This exploration is motivated by the following. The 
He II EW is a measure of the strength of the extreme UV 
(EUV) continuum above 54 eV, compared to the near UV 
continuum, and there are reports that B ALQs have on 
average lower Hen EW th an non-BALQs IjRichards et al.l 
l2002l : iReichard et al.1 [20031 '). A Broad Line Region (BLR) 
wind component is repo rted to be affected by the ioniz- 
ing conti nuum hardness |Leighlv fc Moorell2004l . based on 
Hen EW: lKruczek et al]l201ll . based on aox), and the BLR 



wind component is suggested by I Richards! (|2012f ) to be 
possibly relevant to the BALQ phen omenon. In addition , 
Quv is correlated w ith reddening fe.g.. lBaskin fc Laoiil2005l : 
!Stern fc Laorl!2012! ). and since reddening is more common 
in BALQs (see above), this implies a possible relationship 
between auv and BALQ properties. 

The paper is structured as follows. The data analysis 
method is described in Section [2l In Section [3] we analyse 
composite spectra of the BALQ and non-BALQ samples, 
find a trend between the ionization potential and the av- 
erage BAL depth, and constrain the average H" absorber 
properties. In Section|4]we investigate what parameters span 
the C IV BAL properties. A physical interpretation to our 
findings is proposed in Section [5] In Section [6] we examine 
which dust extinction laws can explain the BALQ reddening 
relative to non-BALQs. Our conclusions are summarized in 
Section [T] 



2 THE DATA ANALYSIS 

The data set is drawn from the SDSS DR7. The object 
BALQ classificatio n is adopted from the !Shen et al.l (l201l[) 
quasar catalogue!^! !Shen et al.l (|201H ) use the [Gibson et al.! 
(2009) BALQ classification for objects that are included in 
the SDSS DR5, and classify the remaining obje cts based on a 
visual inspection of the C iv region. Note that lCibson et al.! 
l!2009r) use a m o dified version of 'balnicity index' (BI) of 
!Wevmann et al.l ([l99l|) to detect BALQs, which they term 
BIq. They integrate the continuum-normalized spectral flux 
starting from a blueshift of 0, rather than —3000 km s~^ 
used in the traditional BI. We include in the data set only 
objects with S/N ^ 3 in the SDSS i-filter, to avoid unusu- 
ally low-S/N spectra (the S/N criterion excludes ~ 5 — 10 
per cent of the objects; see below). We divide the data set 
into two subsets that cover different Arcst ranges as described 
below. 

(i) The 800 ^ Arcst s; 1750 A range i.e., 3.75 < z < 4.25 
for the SDSS (hereafter, the 'high-«' sample). The lower 
limit on z is set to allow a detection of Lyman limit ab- 
sorption intrinsic to the BALQs, and the upper limit is set 
to detect the continuum redward of the C IV emission com- 
plex. The DR7 covers this z range for 228 BALQs and 1320 
non-BALQs. The S/N > 3 criterion leads to 200 BALQs and 
1142 non-BALQs. 

(ii) The 1400 ^ Arcst sS 3000 A range i.e., 1.75 < z < 2.05 
for the SDSS (hereafter, the 'low-z' sample). The lower and 
upper limits are placed to cover the C iv BAL and the 
Mgll emission line, respectively. The DR7 contains 1691 
BALQs and 13,388 non-BALQs in this range. The S/N 
criterion excludes 39 BALQs and 739 non-BALQs. Since 
LoBALQs are a distinct subtype of BAL Qs, with redder 
spectra than the more common HiBALQs (^Wevmann et ajj 
,1991; Spravbcrrv & Fohz 1992; Boroson fc Meyers 1993; 
IReichard et al.!!2003l : IGibson et al.ll2009l ), we exclude from 
the low-2 BALQ sample 56 objects with a detected Mgll 



^ The updated IShen et al.l ll201ll ) quasar catalogue, available at 
|http://d as^sdss.o rg/va/qso_properties-dr7/dr7.htm| includes the 
improved z from lHewett fc Wildl ||201CD . 
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BAL (|Shen et al.ll201lh . and construct a sample of low-z Hi- 
BALQs only. We do not construct a similar high- 2; HiBALQ 
sample, as it is not clear which low-ionization absorption 
line at Are st ^ 1750 A mat ches the Mgll absorption line. 
Note that IShen et al.l (|201ll ) do not conduct a systematic 
search for Mgll BALQs in the post-DR5 quasar sample, and 
report only serendipitously found Mg 11 BALQs for this sam- 
ple. Thus, the exclusion of Lo BALQs from the low -z BALQ 
sample might be incomplete. I Trump et al.l l|2006r ) report a 
Mgll BALQ fraction of 1.3 per cent of quasars for DR3 (i.e. 
164 objects in our sample), where BALs are de tected using 
the 'a bsorption index' AI > criterion, while lAUen et al.l 
l|201ll ). who use the BI > criterion, find a smaller fraction 
of 0.3 per cent for DR6 (38 objects). The different fractions 
result from the different definitions of Al and Bllj Thus, our 
\ow-z HiBALQ sample is likely contaminated by ~ 100 LoB- 
ALQs which pass the AI > criterion, but none of these pass 
the BI > criterion. The final low-z sample is comprised of 
1596 HiBALQs and 12,649 non-BALQs. 

Since the spectra of high-z objects have a relatively low 
S/N (a median value of ~ 7 per resolution element at contin- 
uum regions unaffected by absorption and strong emission jj, 
the spectra are smoothed by a 22 pixel-wide moving aver- 
age filter, which allows to achieve a median S/N of ~ 22 per 
resolution element [the filter width i s equivalent to ^ 9 reso- 
lution elements i.e., ~ 1350 km s"^ iJYork et al.ll2000l ')]. The 
spectra of low-z objects are smoothed by the same filter, 
yielding a median S/N of ~ 66 per resolution element. Since 
we are dealing with BALs (i.e., line width > 2000 km s~^), 
the absorption profiles explored in this study remain well re- 
solved. All spectra are normalized by the mean fiux density 
in the Arest = 1700 - 1720 A rangeQ 

Median spectra are calculated by utilizing a modified 
median method. A simple mean is not adopted as the rep- 
resentative average spectrum, since it can be affected by 
outliers. Outliers have a negligible effect on the median, but 
the median yields a relatively 'noisier' composite compared 
to the mean. We utilize a modified median method, which 
is a hybrid between the standard median and the mean, 
and which mitigates the shortcomings of the two methods. 
First, the object spectra are sorted based on the normal- 
ized /a (Arest). Then, 10 per cent of the objects above and 
10 per cent of the objects below the median are marked. Fi- 
nally, the mean f\ of the marked objects is adopted as the 
composite f\ at a given Arest (a similar method is used to 
calculate the Libor interest rate in Economy). 

The intrinsic BALQ absorption is evaluated by utiliz- 
ing the non-BALQ sample. The non-BALQ composite is 
assumed to represent the intrinsically unabsorbed emission 
of the BALQs (the construction of composites matched in 



^ AI > (BI > 0) requires at least one continuous absorption 
trough with a minimal absorption depth of 0.1 and a minimal 
width of 1000 km s'^ (2000 km s'^) in the -29,000 < I'shift < 
km s^i (-25, 000 < -Ushift < -3000 km s"!) range. 
^ The S/N of an object is estimated by calculating the ratio 
between the mean and standard deviation of fx in the Arest = 
1700 - 1720 A window. 

■* It should be noted that the wavelength used for the normaliza- 
tion is not significant, as verified by repeating the main analysis 
using the Arest = 1275 — 1285 A range instead. 



emission properties is described below). The ratio between 
the BALQ and non-BALQ composites (hereafter the R spec- 
trum) is interpreted as the intrinsic BALQ absorption spec- 
trum. Since the ratio does not remove an overall spectral 
energy distribution (SED) differences between BALQs and 
non-BALQs, we fit a straight-line continuum across the pro- 
file from Ushift = to —30, 000 km s"'^. A residual Siiv emis- 
sion, if present, has little effect on the C iv measurements. 
The caveat of our approach is that if there are significant 
intrinsic differences in line emission within the absorbed re- 
gion for BALQs, then the derived absorption profile will be 
biased. 

The individual absorption profiles are often character- 
ized by a relatively narrow absorption peak. This charac- 
teristic profile is smeared out when forming the composite, 
due to the distribution of peak absorption velocities. In or- 
der to investigate the median 'BAL rest-frame' absorption 
profile near the absorption peak, we also form a composite 
based on the aligned C iv peak-absorption. This composite, 
calculated for the high-z BALQs, allows to detect relatively 
narrow (FWHM > 1500 km s~^) and weak absorption lines 
which are otherwise heavily blended in the crowded FUV re- 
gion, and also provides some hint on the optical depth of re- 
solved multiplet absorption. The BALQ spectra are aligned 
by shifting each spectrum in velocity, so that the maximum 
absorption of C iv falls at the same velocity for all spec- 
tra. Then, a composite of the aligned spectra is calculated. 
The 22-pixel spectral smoothing by a relatively broad filter, 
smooths out any narrow features (< 500 km s^^) super- 
imposed on the C iv BAL, and reduces their effect on the 
alignment procedure. When the matched non-BALQ sam- 
ple is constructed, a similar distribution of velocity shifts is 
used. 

We study the dependence of the BALQ absorption 
properties on the following parameters. 

(i) The Hell EW, which may serve as a measure of the 
EUV ionizing SED hardness. The Hell EW is measured by 
integrating the normalized fx in the Arest ~ 1620 — 1650 A 
range, assuming fx"^^ ~ 1 i-e., a constant /a°"' between 
the normalization window (1700-1720 A) and ~ 1620 A. 
The adopted Arest range minimizes contributions from C iv 
and Oiii] AA1661, 1666 to the evaluated Hen EW. Using 
a constant /a°"* corresponds to a spectral slope of a = —2 
(/y oc J/"), which serves as a rough approximation for the lo- 
cal slope (Section l4.2|l . A significantly more accurate deriva- 
tion of the local /™"* is hindered by the Civ BAL for 
BALQs. The derived Hell EW is slightly overestimated in 
objects where the continuum slope is bluer, but these objects 
tend to have a higher He 11 EW, where the exact continuum 
placement is less significant. The constant />°"' is a good 
approximation in the low Hell EW objects. In the reddest 
objects, the continuum at Arest < 1700 A drops below the 
adopted constant /x°"', leading to slightly negative He 11 EW 
(see Section r4.3|l . The small errors in the absolute values of 
the Hell EW are significantly smaller than the trends with 
the Hell EW explored below. 

(ii) The UV spectral slope auv, which may serve as a 
measure of dust absorption. For the high-z sample, the slope 
is evaluated between the 1275-1285 and 1700-1720 A win- 
dows, and is denoted as auvs. The latter window is redward 
of Ollll A1665 emission and the Civ emission 'shelf, and 
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blueward of N ill] A1750. The former window is between the 
red wing of Ly a+N v emission complex and O I A1303. The 
window between Ol A1303 and Cll A1335 is not utilized, 
because it is affected by C ll and Si iv absorption in BALQs. 
For the low-z sample, the slope is evaluated between the 
1700-1720 and 2990-3010 A windows, and is denoted as 
Quvi- The latter window is redward of Mgll A2798. 

We also explore the dependence of the BALQ absorption, for 
the low-2 sample, on the following additional parameters. 



(i) L(3000A). taken from lshen et al.1 (l201ll). 

(ii) Mgii FWHM, taken from[Shene^3'^oi3). 

( iii) M bh , estimated using parameters (i) and (ii), and 
the Vestergaard fc Osmer (2009 ) prescription (values listed 
in lShenet al.ll201ll ). 

(iv) L/L^dd, evaluated as logL/L^dd = logL(3000A) — 
log Mbh — 3 7.4, where we ad opt a bolometric correction fac- 
tor of 5.15 (|Shen et al.ll2008l ). 

We do not use the C iv emission line to estimate A^bh given 
the significant uncer tainties associated with this line (e.g. 
iBaskin fc Laorll2005l . and citations thereafter). The objects 
are binned for each parameter, so that each BALQ bin con- 
tains the same number of objects. The corresponding non- 
BALQ bins span the same parameter range, but the distri- 
bution of objects within each bin may be different. A median 
composite spectrum is calculated for each bin. 



3 THE AVERAGE BALQ PROPERTIES 

Figure [l] presents a comparison between the high- 2; BALQ 
and non-BALQ composites, aligned by their tabulated z. 
We present both the median and the mean composites, to 
demonstrate the level of systematics produced by the av- 
eraging procedure, which is larger than the purely statisti- 
cal noise. The two composites are similar overall, except at 
wavelength regions of strong absorption and emission fea- 
tures (e.g., Lya), where the mean composite is more af- 
fected by outliers. Bottom panel presents the R composite. 
The median and mean R composites are similar. There are 
small differences in absorption amplitude, which are likely 
selection effectslj We indicate the laboratory wavelength lo- 
cation of absorption lines of various ionized species. There 
are many absorption lines in the Arost ~ 800 — 900 A range, 
and the indicated absorption line locations are for illustra- 
tion only. The steep decline of flux blueward of the Lyman 
edge for both BALQ and non-BALQ composites (top panel) 
is caused by intervening Lyman a systems. Observations 
of low-z AGN, which are unaffected by i ntervening Lyman 
systems, do not show this steep decline (IZheng et al.l 1 19971 : 
iTelfer et al.ll2002l : IScott et al.|[20oi : IShull et al.ll2012l ). 

The three main features which can be seen in the R 
composites are: (i) there is no detectable Lyman limit ab- 
sorption; (ii) the high ionization lines observed shortward of 
Lya, in particular Nv and Ovi, show stronger absorption 



^ The reversal of the mean versus median _R spectra at the C IV 
and Ovi absorption troughs is not a significant result. It reflects 
the fact that the residual fiux at the O VI trough is closer to zero 
than for C IV. 



than observed longward of Lya; (iii) the overall continuum 
ratio is red. These results are further discussed below. 

Figure [5] compares the absorption profiles of Siiv, Civ, 
N V and O vi. The higher the ionization potential, the 
stronger is the BAL trough. The Silv absorption line, pro- 
duced by photons above 33.5 eV and destroyed by photons 
above 45.1 eV, has the shallowest trough, ft is followed in 
trough depth and absorption EW by Civ, Nv and Ovi 
(47.9-64.5, 77.5-97.9 and 113.9-138.1 eV, respectively). The 
absorption EW (transmission at maximal absorption) is 5.8 
(0.87), 9.3 (0.83), 24.6 (0.56) and 25.2 A (0.46) for Siiv, 
Civ, Nv and Ovi, respectively. The Nv absorption may be 
affected by Lya at tiahift = 0. Note that Cll A1335 falls on 
the same trend, as its absorption is weaker (EW ^ 1 A) 
than Si IV (Fig. [U bottom panel). 

Our findings are consistent with the I Allen et al.l (|201ll ) 
report that on averag e Siiv reaches small er val ues of peak 
absorption than C iv. iGibson et al.l (|2009|) and lAllen et al.l 
(2011) find that the Civ BI distribution has a marginally 
larger tendency towards higher values than Si iv. The differ- 
ence in absorption depth likely implies difference in CF. 

Figures [3] and [4] present the C iv peak-absorption 
aligned R spectrum and absorption line profiles, respec- 
tively. Line identifications are indicated in Fig. [3] at the 
line laboratory wavelength. The peak-absorption alignment 
method aligns all the BALs, including low-ionization BALs 
(e.g., Cll A1335). The line profiles of the BALs are exam- 
ined further in Fig. [l] where the lines are grouped vertically 
based on the ionization potential. The absorption profiles are 
dominated by the relatively narrow absorption component 
(FWHM ^ 3000 km s"^). Fig. H indicates that for several 
lines the absorption is saturated, as indicated by the similar 
absorption depth of doublet components with different oscil- 
lator strengths. The absorption profile in these cases is set by 
a velocity dependent CF (e.g., S vi AA933, 945). Other lines 
are not saturated (e.g., Siiv AA1397, 140 3), as can be indi- 
cated by more quantitative analysis (e.g. iDunn et al.ll2012l . 
Capellupo et al. in prep, and Hamann et al. in prep.). 

3.1 Constraints on the BAL H" Column 

Figure 3] can be utilized to constrain the average BAL H" 
column A/"(H°), based on the Lya absorption line. We begin 
by deriving A'^(H'') assuming a uniform foreground screen 
i.e., CF = 1. The Lya absorption line resides within the 
Nv absorption trough. We estimate the Lya absorption 
depth by interpolating linearly between the two 'shoulders' 
at «shitt = ±2000 km s^^ (Fig. |3|, which yields a contin- 
uum fx ~ 0.55 at the line centre. The measured fx ~ 0.42 
at the peak absorption implies a Lya peak absorption of 
exp(— To) = 0.42/0.55 ~ 0.76. Assuming a Maxwellian ab- 
sorption profile with a velocity dispersion b para meter of 
1000 km s"\ yields iV(H'') ~ 1 x 10^^ cm"^ ('e.g.. lOrainj 
l201li . eq. C.l). This N(R°) implies r ~ 0.01 at the Lyman 
edge (e.g. JOsterbrockll 19891 . eq. 2.4), which is consistent with 
a non-detection of the edge in the absorption-peak aligned 
R spectrum (Fig. 3] Ly lim. panel). However, the Lya line 
and Lyman edge can be highly saturated, with A^(H'') well 
above the derived value. In that case, the Lya absorption 
profile reflects the CF of the saturated cold absorber, which 
is at least 1 — 0.76 = 0.24. However, the observed Lyman 
edge region shows no detectable edge, i.e. r < 0.1, and im- 
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Figure 1. Comparison between the BALQ and non-BALQ composite spectra. Top panel: the median normalized flux of the BALQ 
and non-BALQ samples (thick and thin black solid line, respectively). The mean normalized flux is also presented for comparison (red 
lines). The composites can differ at strong absorption and emission features (e.g., Lya), where the mean is more affected by outliers. 
The Lyman limit (vertical dotted line) and prominent emission lines are indicated. Bottom panel: the ratio between the BALQ and 
non-BALQ composites. The laboratory wavelength location of prominent absorption lines is indicated longward of the Lyman limit. We 
also indicate for reference the location of possible lines below the Lyman limit. The indicated Sll* and Sll** lines have excitation energy 
of 1.8 and 3.1 eV, respectively. There is no significant absorption edge detected at the Lyman limit. Note also the declining continuum 
ratio to the blue, and the increasing absorption features depths. 



plies iV(H'') < 1 X 10" cm"^ if CF = 0.24. In addition, 
the Ly /3 shows only weak absorption, well below 0.24, if we 
adopt the O vi shoulder as the local continuum, which also 
suggests a non saturated Lya absorption. 

Another constraint on the CF of any nearly neutral sat- 
urated absorber can be derived from the Cll A1335 absorp- 
tion line. The R spectrum in the C ll region shows f\ ~ 0.84 
at the line centre, compared to an estimated local continuum 
of 0.89 (Fig.Bl Cll 1336 panel), which implies a CF < 0.06, 
consistent with the non detection of the expected Lyman 
edge absorption. 

We note in passing that photoionization models indicate 
that the Ly a absorption should be at least as strong as Si iv 
for solar metallicity, regardless of the ionization parameter 
(jBaskin &: Laoill2012l ). This is consistent with our observa- 
tions, as the Si IV absorption goes down from f\ « 0.9 to 



0.62 that indicates peak absorption of ~ 0.3, which is simi- 
lar to the Lya peak absorption. 

An obvious general caveat is that the above analysis 
is applied to the median absorption profile, and may not be 
valid to individual objects, such as LoBALQs, which are not 
excluded from the high-z BALQ sample. 



4 WHICH PARAMETERS CONTROL THE 
ABSORPTION PROFILE? 

In this section we use composite spectra to examine the 
relationship of Civ BAL properties to Hell EW, auv, and 
other measured parameters. 
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Figure 2. The absorption profiles of the prominent absorption 
lines SiV, Civ, Nv and Ovi. The velocity scale is set by the 
longer wavelength of each doublet. The profiles are normalized to 
1 at Vshiit = km s~^. Laboratory wavelength location of various 
other possible absorption lines is indicated by vertical tick marks 
(with the same colour coding). The absorption depth increases 
with the ionization potential of the absorbing ion. 



4.1 Civ BAL trends with the Hen EW 

Figure [5] presents the liigii-z BALQ and non-BALQ samples 
separated into two bins based on their He ll EW (as furtiier 
described below). The upper panel presents the composite 
spectra of the two bins for BALQs and for non-BALQs, and 
the lower panel presents the R spectra for the two bins. The 
low Hell EW bin displays a somewhat deeper and signifi- 
cantly broader C iv BAL, with peak absorption shifted to 
larger Ushift- This trend is also observed for Siiv and possi- 
bly for O VI, although a clear detection of the trend for the 
latter is hindered by the low S/N and the blending of a num- 
ber of adjacent absorption lines. The trend does not appear 
significant for N v, probably due to contamination by the ad- 
jacent low ionization line Lya, which is prominent relative 
to other low ionization lines (e.g., Fig.[3|. The weaker trend 
between the He ll EW and the N v and O vi BALs is further 
discussed in Section [5. II The drop in the R spectrum at the 
red wing of the C iv emission line reflects a mismatch in the 
C IV emission EW between BALQs and non-BALQs for the 
same bin (see Fig. [S] top panel) , and is most likely not an 
absorption effect. Table [T] summarizes various properties of 
the two He ll EW bins, and lists the first three moments of 
the C IV absorption profile (EW, mean Wshift and dispersion 
a). 

Figure [6] presents the He ll EW binning for the low- 
z sample, this time binned to 4 bins (as further described 
below), allowed by the larger sample size (see Section ^. 
The same trend of increasing depth and blueshift in the 
C IV absorption profile with decreasing He ll EW is ob- 
served here as well. There is also an excess emission in the 
C III] -I- Si III] -f Al III complex in the 7? spectra, which does 
not depend on the He ll EW. The excess is related to a trend 
with auvi discussed below (Section 14. 2[) . Table [5] lists the 
bin median properties, including the median Mgll FWHM, 
I/(3000A), Mbh and L/Z/Edd- It also presents the number of 



LoBALQs with similar Hell EW. LoBALQs are found pre- 
dominately in the lowest He ll EW bin (43 out of 56 objects). 
In Appendix \K\ we present a LoBALQ composite based on 
the 56 objects we excluded from our low-z sample in order 
to form the low-z HiBALQ sample (Section [2|. The redder 
R spectrum of the lowest Hell EW bin, compared to the 
other three, might represent a smooth transition between 
HiBALQs and LoBALQs, where the latter have the low- 
est Hen EW and the reddest R spectrum (Appendix E| . 
The smooth transition scenario can be further explored by 
analysing the Mgll BAL dependence on the Hell EW and 
reddening in LoBALQs. This requires a larger sample of 
LoBALQs and is outside the scope of this paper. 

Tables [T] and [2] list for each BALQ bin the observed 
BALQ fraction in DR7, which is defined as /balQs = 
A'^balQs/(A''balQs + A'non-BALQs)- The tabulated errors on 
/balQs are minimal, and are based purely on number statis- 
tics. There is a strong trend between /balQs and the Hell 
EW. For the low- 2; sample, /balQs goes up from 7.0 ± 0.4 
per cent for the highest Hell EW bin, to 23 ± 1 per cent 
for the lowest bin. For the high- 2: sample, a similar trend is 
found, and /balQs goes up from 11 ± 1 per cent for the high 
Hen EW bin to 22 ± 2 per cent for the low Hell EW bin. 



4.2 C IV BAL trends with auv 

Figures[7]and[S]explore the dependence of the BALQ absorp- 
tion on auv- Figure[7]presents the high-2 sample binned into 
two auva bins. The Civ absorption profile follows the trend 
observed in the He 11 EW binning, where the redder auvs bin 
corresponds to the lower Hell EW bin. Table [3] lists the bin 
median properties of the high-2 sample. The median auvs 
of the BALQ bins is redder than the slope of the matched 
non-BALQ bins. Although the bins are matched by the auvs 
range (Section[2|, the distribution of auvs values within the 
bin is different for the BALQs and non-BALQs. This slope 
difference yields the observed residual reddening of the R 
spectrum (Fig. [3 bottom panel). Figure [S] presents the low- 
z sample divided into four Quvi bins. The trend with auvi 
is different from the trend with auvs- The low velocity Civ 
absorption profile becomes deeper as auvi becomes redder, 
while at higher velocities the C iv absorption profile remains 
unchanged. Table |4] summarizes the bin median properties. 
Note that the C ill] -I- Si ill] -I- Al ill complex emission becomes 
stronger with redder auvi (Fig.|Hl top panel), and the BALQ 
and non-BALQ matched composites have a similar complex 
emission. 

The fraction of BALQs increases as the objects get red- 
der. At the high-z sample, /balQs increases from 11 ± 1 to 
26 ± 3 per cent, while at the low-2 sample, /balQs increases 
from 6.9 ± 0.4 to 21 ± 1 per cent from the bluest to the red- 
dest quartile. Most of the detected LoBALQs (44 out of 56) 
have slopes at least as red as the reddest bin, as expected 
(see Section [T|. 



4.3 Civ BAL trends with both Hell EW and auv 

Figure [9] shows a zoom-in on the C iv absorption profile 
for the low-z sample, and highlights the strong and dif- 
ferent trends between the C iv BAL properties and the 
Hen EW and auvi (not true for auvs)- The Hell EW 
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Figure 3. The absorption profiles near tlie absorption peak. The spectra are aUgned so that the C IV absorption dip falls at 1549.48 A 
(i.e., a 1:1 contribution to absorption from the two components of Civ AA1548.20, 1550.77). Note that both the low- and the high- 
ionization BALs are at the same Ughift a-s Civ. This allows to separate heavily blended features, such as Lya and Nv, and to detect 
weak and relatively narrow features, which are otherwise heavily blended, such as Svi (Fig.[lJ. 



Table 1. 


The median properties 


of the He 


II EW binned hig 


h-z objects. 




Class 


Bin 


Nobj 


Hen 

EW 

[A] 


auvP 


EW 

[A] 


C IV BAl!^ 

[km s-i] [km s-i] 


/balqJSI 
[per cent] 


BALQs 


high 

low 


100 
100 


5.4 
0.4 


-1.04 
-1.76 


12.0 

14.8 


-13000 
-12500 


6800 
4700 


11 ±1 
22 ±2 


non- 
BALQs 


high 

low 


778 
353 


6.1 

0.7 


-0.71 
-1.36 











" Measured between 1280 and 1710 A {f„ oc u"). 

Measured between iighift = and —30,000 km s~^. 
"^ Measured from the total quasar population. Errors are based on number statistics. 



Table 2. The median properties of Hen EW binned low-2 objects. 



JT 



/balqP 
[per cent] 



Class 



Nbin 



N, 



obj 



Hen 

EW 

[A] 



auv 



P 



L(3000A) 
[erg s-i] 



Mgll 
FWHM 

[km s"-"-] 



A/bh 

[Mq] 



log 
L/L-Edd 



A^LoE 



C IV BAlJSl 



EW 

[A] 



shift 

[km s^-"-] [km s^-"-] 



BALQs 



399 
399 
399 
399 



7.0 
4.7 
3.0 
0.9 



-0.97 
-0.83 
-0.79 
-0.96 



45.90 
45.97 
46.00 
45.99 



3900 
4400 
4100 
3400 



9.03 
9.14 
9.09 
8.93 



-0.48 
-0.52 
-0.44 
-0.30 



7 
1 
5 
43 



7.8 
13.9 
15.7 
22.5 



-4600 
-9200 
-10900 
-11900 



2500 
5700 
5900 
6000 



7.0 ±0.4 

10.3 ±0.5 

14.1 ±0.8 

23 ±1 



non- 


1 


5330 


7.3 


-0.77 


45.81 


3800 


8.93 


-0.49 


BALQs 


2 


3461 


4.7 


-0.65 


45.93 


4000 


9.03 


-0.47 




3 


2425 


3.1 


-0.61 


45.94 


3700 


8.98 


-0.38 




4 


1368 


1.0 


-0.69 


45.93 


3300 


8.87 


-0.30 



" Measured between 1710 and 3000 A{fi, ocu"). 

Number of LoBALQs (identified based on a detection of Mgll BAL: IShen et al.ll2011l ) with Hell EW within or nearest to the bin range. 
"^ Measured between tighift = and —30,000 km s~^. 
" Measured from the total quasar population. Errors are based on number statistics. 
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Figure 4. Line profiles of tlic aligned composite spectrum (see Fig.[3]|. The lines are grouped according to the ionization potential, and 
are arranged from the longer to the shorter wavelength at each group. The velocity scale is calculated relative to the redder line of each 
multiplet, which is indicated at each panel. A zero velocity shift and fghift ol other absorption lines of each multiplet are marked (dashed 
and dotted line, respectively). The location of other absorption lines is indicated by vertical tick marks at the top of each panel. Note 
that some of the doublets appear optically thick (e.g., S VI 945), while others appear unsaturated (e.g.. Si IV 1403). 



Table 3. 


The median 


aroperties 


of auvs 


binned high-2 ob 


ects.lSl 




Class 


Bin 


A^obj 


"UVs 


Hen 

EW 

[A] 


EW 

[A] 


CivBA 

.,mcan 
^shift 

[km s'l] 


[km s-i] 


/balqP 
[per cent] 


BALQs 


blue 
red 


100 
100 


-0.92 
-1.90 


4.8 
1.2 


7.9 

17.2 


-11700 
-12600 


6300 
4800 


11 ±1 
26 ±3 


non- 
BALQs 


blue 
red 


848 
282 


-0.72 
-1.66 


5.3 

1.2 











"• Measured between 1280 and 1710 A (/,. oc i/"). 

" Measured between i>shift = and —30,000 km s~^ . 

^ Measured from the total quasar population. Errors are based on number statistics. 
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Figure 5. The dependence of the intrinsic absorption on the Hell EW for the high-^ sample. Top panel: the composite spectra of 
high and low Hell EW BALQs and non-BALQs. The low Hell EW spectra are shifted down by 0.7 for presentation purposes. Bottom 
panel: the ratio between the corresponding BALQ and non-BALQ bins. The ratio of the complete samples is also presented (same as in 
Fig-EJi which shows an apparent red side absorption in Civ that results from the different median Civ emission profiles in BALQs and 
non-BALQs. This effect is partly compensated by using the two subsamples, which have better matched C IV emission strengths. Note 
the significantly stronger C IV absorption for the lower He II EW composite. A similar effect in seen in the S IV BAL, some effect may be 
seen in the Ovi BAL, and only a marginal effect is seen in Nv, possibly due to blending with the low ionization Lyo BAL. 



mostly controls the characteristic absorption Ushift, while 
auvi mostly controls the absorption depth, in particular at 
Vshift < 10,000 km s~^. Thus, it appears that the Hen EW 
and auvi are independent parameters which span the C iv 
absorption profile properties. Their independence can also 
be inferred from Table [21 which shows similar median auvi 
values in the different Hell EW bins, and similarly in Ta- 
ble m where the different auvi bins show similar median 
Hell EW values. This is in contrast with the high-z sample, 
where auvs and the He ll EW are correlated (Tables [T] and 

EJ. 

The highest Hen EW and the bluest auvi bins may be 
incomplete in other BALQ studies based on a BI selection. 
A value of BI > requires a BAL with a significant absorp- 
tion (> 10 per cent) at least down to lighift = —5000 km s^^ 
(given a minim al width of 2000 km s~ ^, starting at Ushift = 
-3000 km s~\ IWevmann et al.|[l99lf) . Since we find that 



t^shift decreases with increasing Hen EW, the highest Hen 
EW bin might be missing yet lower fahift BALQs, with peak 
absorption at i^shift > —3000 km s~^. Similarly, BALQs re- 
siding in the bluest auvi bin may be missed, since the peak- 
absorption is < 10 per cent for this bin (Fig. [9]). Our sample 
is derived based on a Bio and visual selection (see Section[21), 
where Bio > requires significant absorption (> 10 per 
cent) at least down to Wahift = —2000 km s~^. Thus, large 
He II EW, low tfshift BALQs are probably not missing in our 
sample, but weak BALs (< 5 per cent) possibly present in 
BALQs with the bluest auvi niay be missing. 

Figure [TO] explores the C iv BAL dependence on the 
Hen EW for different values of auvi- We divide the \ow-z 
sample into 4x4 bins based on auvi and the He n EW. First, 
the BALQ sample is divided into four auvi 'parent' bins, 
and matched non-BALQ bins are constructed (same proce- 
dure as in Fig. [8]). Then, each auvi parent bin is divided 
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Figure 6. Same as Fig. [3] for the low-z sample. The BALQ and non-BALQ samples are divided here to four subsamples. The 2nd to 
4th bin spectra in the upper panel are shifted down for presentation purposes by 0.6, 1 and 1.3, respectively. The Civ absorption profile 
becomes stronger and broader with lower Hell EW values, as seen in the high- 2; sample (Figure [SJ. 



Table 4. The median properties of auvi binned low-2 objects.!^ 



Class 



NhlYL Wobj CtUVl 



Hen 
EW 

[A] 



L(3000A) 



Mgii 
FWHM 

[km s^^] 



log 
Mbh 

[M(d] 



log 



Nl^a 



JT 



EW 

[A] 



C IV BAI0 

shift 

[km s^-"-] [km s" 



/balqP 

[per cent] 



BALQs 


1 


399 


-0.48 


3.5 


45.92 


4100 


9.08 


-0.49 


3 


7.6 


-10900 


6200 


6.9 ±0.4 




2 


399 


-0.76 


4.0 


45.97 


4100 


9.09 


-0.45 


3 


8.2 


-9000 


5600 


10.7 ±0.6 




3 


399 


-1.05 


4.3 


45.97 


3900 


9.05 


-0.41 


6 


6.8 


-4700 


2300 


14.5 ±0.8 




4 


399 


-1.51 


3.5 


46.00 


3600 


9.02 


-0.38 


44 


10.3 


-4500 


2500 


21 ±1 


non- 


1 


5414 


-0.44 


4.6 


45.86 


3800 


8.96 


-0.46 












BALQs 


2 


3334 


-0.76 


5.4 


45.92 


3800 


8.99 


-0.43 














3 


2349 


-1.02 


5.9 


45.89 


3700 


8.97 


-0.42 














4 


1526 


-1.45 


5.1 


45.86 


3600 


8.92 


-0.41 













" Measured between 1710 and 3000 A (f^ oci/°'). 

Number of LoBALQs (identified based on a detection of Mgll BAL: IShen et al.ll201lh with auvi within or nearest to the bin range. 
^ Measured between tighift = a-nd —30,000 km s~^. 
" Measured from the total quasar population. Errors are based on number statistics. 
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Figure 7. The dependence of the intrinsic absorption on the «uvs for the high-z sample. Top panel: the composite spectrum of BALQ 
and non-BALQ subsamples matched in ouva range. The BALQ and non-BALQ red bins are shifted down by 0.5 for presentation 
purposes. Bottom panel: the ratio between the corresponding BALQ and non-BALQ composites. Note the stronger absorption for the 
redder QuVs- There is a residual reddening of the R spectrum, since the median slope of the BALQ bin is redder than the slope of 
matched non-BALQ bin. 



into four Hen EW bins, with matching non-BALQs from 
the corresponding auvi non-BALQ bin. As Fig. [10] shows, 
the Civ absorption at a given ayvi, shifts to higher fahift 
as the Hell EW decreases, at all four auvi bins. Also, the 
maximum absorption depth at all He II EW bins, gets larger 
as Quvi gets redder. The change in Vshift is most prominent 
for the bluest two bins. Table [5] lists the median properties 
for the 4x4 bins. Note that the reddest and lowest Hell 
EW BALQ and non-BALQ bins have median He ll EW < 0. 
This is an artefact which results from measuring the EW by 
using the continuum at 1700-1720 A, where in the reddest 
and weakest Hell objects the flux density at Hell is lower 
than at 1700-1720 A, leading to a negative EW (Section [2]). 
The negative Hell EW has no effect on the binnings, which 
are independent of the absolute He ll EW values. The value 
of /balQs increases with decreasing Hell EW, at a given 
Quvi, and also as Quvi gets redder, at a given Hell EW. We 
find the smallest value of /balQs ~ 4.2 ± 0.4 per cent in the 
highest Hell EW and bluest Ofuvi bin, increasing to 31 ± 4 



per cent in the lowest He ll EW and reddest auvi bin. Most 
of the LoBALQs (33/56) reside in the reddest and lowest 
Hen EW bin. 



4.4 C IV BAL dependence on other emission 
properties 

Figure [TT] investigates whether a certain physical property 
underlies the above relations between the absorption profile 
and the Hell EW and Quvi. The figure presents the Civ 
BAL profile for samples binned based on the Mgll FWHM, 
I/(3000A), A/bh and L/L-^dd values. Table [6] lists the prop- 
erty median value of each bin for each of the 4 binning pro- 
cedures. There is no trend in the absorption profile with 
I/(3000A). There is a trend of stronger and broader absorp- 
tion profile with decreasing Mgll FWHM, however the dy- 
namical range in the absorption properties is significantly 
smaller than the range found for the He ll EW or auvi bin- 
ning. Similar trends are seen with Mbh and L/L-Edd- How- 
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Figure 8. The dependence of the intrinsic absorption on «uvi foi' the low-z sample. The BALQ and non-BALQ samples are divided 
into four bins. The 2nd to 4th bin spectra are shifted down for presentation purposes by 0.6, 1.3 and 2, respectively. The C IV absorption 
profile becomes deeper as «uvi becomes redder at the lower velocities, and remains unaffected at the highest velocities. 



ever, these two parameters are strongly correlated with the 
Mgll FWHM, given the small range in L(3000A) (Table©, 
and since I/(3000A) displays no trend, the C iv BAL trend 
with Mbh and L/I/Edd is equivalent to the trend with Mg ll 
FWHM. Also, /balq is either nearly constant, or shows a 
range of generally less than a factor of two, which is another 
indication these four parameters are not the fundamental 
parameters that underlie the C iv BAL profile. Thus, the 
Hen EW and auvi appear to be the primary parameters. 
The possible underlying physical explanation is discussed 
below (Section [5]). 



A cautionary note is that one should keep in mind that 
the evaluated trends with Mbh and L/L^aa may not appear 
significant, since the sample does not cover a large enough 
dynamical range in L (> 1.5 dex) to surpass the ~±0.5dex 
uncertainty i n the Mbh prescr i ption calibration (JLaoill 19981 : 
lKrolikll200ll : iBentz et al.ll2009l : Iwoo erai]|2010l ). 



4.5 Comparison with previous studies 

The resu lts of this study a re con sistent with previous find- 
ings. Ri chards et all (|2002l ) and iReichard etaP ()2003l ) re- 
port that B ALQs have sligh tl y wea ker Hen emission than 
non-BALQs. iRichards et all |2002l ) find that their LoB- 
ALQ composite has a weaker Hen emission than the Hi- 
BALQ and non-BALQ composites. This is consistent with 
most of the Lo B ALQ s residing in our lowest Hen EW bin. 
iReichard et all (J2003l ) also find the HiBALQs to be red- 
der than non-BALQs , and the LoBALQs to be the reddest. 
[Richards et al.l l|2003h report a trend between the quasar 
colour and /balQs by similarly dividing their sample into 
four bins. Their /balQs goes up from ~ 3 per cent for the 
bluest bin (compared to our 6.9 ± 0.4 per cent) to ~ 18 
per cent for the reddest bin (21 ± 1 per cent). A relation- 
ship between reddening and the Hen EW that breaks for 
the reddest bin (Tab le |4]) was also previously reported by 
[Richards et al.l (|2003t ). The large fraction of LoBALQs found 
here in b ins that have the l argest C IV absorption is consis- 
tent with lAllen et all (|201ll ). who find that LoBALQs reside 
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Table 5. The median properties of auvi a-nd the Hell EW binned low-2 objectslSl 

A^LoBAlJ^ 



/balqJ^ 
[per cent] 



Class/ 
"UVl 



A'bin A^obj 



Hen 
EW 

[A] 



log 
L(3000A) 

[erg s-i] 



Mgll 
FWHM 

[km s-i] 



log 



log 

L/L-Edd 



C IV BAlP 



EW 

[A] 



shift 

[km s^-"-] [km s"" 



BALQs 


1 


99 


6.1 


45.84 


4500 


9.08 


-0.63 


2 


6.2 


-6000 


4200 


4.2 ±0.4 


-0.48 


2 


100 


4.2 


45.91 


4700 


9.17 


-0.61 





12.2 


-10200 


6000 


6.3 ±0.6 




3 


100 


2.9 


45.98 


3900 


9.06 


-0.42 





13.9 


-13000 


6300 


9.0 ±0.9 




4 


100 


1.2 


45.94 


3300 


8.90 


-0.32 


1 


19.5 


-13500 


5600 


14 ±2 


non- 


1 


2276 


6.5 


45.77 


3900 


8.95 


-0.54 












BALQs 


2 


1488 


4.3 


45.91 


4000 


9.03 


-0.47 












-0.44 


3 


1015 


2.9 


45.93 


3600 


8.98 


-0.37 














4 


598 


1.2 


45.90 


3200 


8.84 


-0.30 












BALQs 


1 


99 


6.5 


45.90 


4600 


9.14 


-0.59 





6.3 


-4200 


1800 


6.1 ±0.6 


-0.76 


2 


100 


4.7 


45.97 


4300 


9.12 


-0.51 





13.0 


-10500 


6200 


11 ±1 




3 


100 


3.3 


46.04 


4300 


9.16 


-0.50 


1 


15.9 


-10400 


5300 


15 ±2 




4 


100 


1.3 


45.97 


3400 


8.93 


-0.29 


2 


23.7 


-13700 


6200 


24 ±3 


non- 


1 


1518 


6.9 


45.88 


3800 


8.98 


-0.46 












BALQs 


2 


827 


4.7 


45.96 


4000 


9.05 


-0.44 












-0.76 


3 


547 


3.4 


45.97 


3700 


8.99 


-0.38 














4 


318 


1.6 


45.97 


3400 


8.92 


-0.27 












BALQs 


1 


99 


7.7 


45.91 


3700 


8.95 


-0.43 





8.5 


-4200 


2000 


8.2 ±0.9 


-1.05 


2 


100 


5.1 


45.97 


4300 


9.18 


-0.51 





11.8 


-6200 


3400 


16 ±2 




3 


100 


3.4 


46.02 


3900 


9.08 


-0.41 


1 


14.2 


-8800 


5400 


19 ±2 




4 


9P 


1.2 


45.97 


3400 


8.95 


-0.31 


5 


22.1 


-10900 


5900 


29 ±3 


non- 


1 


1107 


8.0 


45.83 


3700 


8.93 


-0.44 












BALQs 


2 


539 


5.1 


45.95 


3900 


9.04 


-0.42 












-1.02 


3 


414 


3.6 


45.95 


3800 


9.00 


-0.41 














4 


243 


1.3 


45.97 


3600 


8.95 


-0.31 












BALQs 


1 


99 


7.7 


45.93 


3300 


8.86 


-0.37 


3 


6.3 


-3100 


1400 


14 ±1 


-1.51 


2 


100 


4.7 


46.05 


4000 


9.09 


-0.46 


3 


18.3 


-6900 


4700 


22 ±2 




3 


100 


2.2 


46.02 


3700 


9.05 


-0.37 


5 


14.6 


-5700 


3700 


25 ±3 




4 


100 


-0.6 


46.03 


3700 


8.93 


-0.33 


33 


25.9 


-9600 


5800 


31 ±4 


non- 


1 


620 


8.4 


45.78 


3500 


8.85 


-0.44 












BALQs 


2 


362 


4.8 


45.91 


3900 


9.01 


-0.45 












-1.45 


3 


307 


2.4 


45.95 


3500 


8.93 


-0.37 














4 


224 


-0.2 


45.91 


3400 


8.87 


-0.31 













" The slope is measured between 1710 and 3000 A (/^ oc 1'°'). 

^ Number of LoBALQs (identified based on a detection of Mg II BAL; IShen et al.ll2011f) with auvi and the He II EW within or nearest to 

the bin range. 

'' Measured between Dghift = and —30,000 km s~^. 

Measured using the total quasar population. Errors are based on number statistics. 
^ An additional BALQ has Hell EW smaller by more than 1 A than the minimal Hell EW of the corresponding non-BALQ bin, and is 
excluded from the BALQ bin. 

Table 6. The property median value of Mgll FWHM, L(3000A), Mbh and L/L^dd binned low- 2; samples^ 



Class 



A^h 



[km s" 



Mgll FWHM 

'^] -^obj /bALQs 



[erg s- 



logL(3000A) 
-'] A^obj /bALQs [Mq] 



log Mbh 

^obj /bALQb 



log L/L-Edd 

-^obj /bALQb 



BALQs 



6200 
4500 
3400 
2200 



399 
399 
399 
399 



13.5 ±0.7 
11.2 ±0.6 
9.8 ±0.5 
11.1 ±0.6 



46.26 
46.04 
45.90 
45.66 



399 
399 
399 
399 



15.2 ±0.8 
13.9 ±0.7 
11.9 ±0.6 

7.4 ±0.4 



9.47 
9.16 
8.91 
8.49 



399 
399 
399 
399 



15.5 ±0.8 
12.2 ±0.6 
9.4 ±0.5 
9.7 ±0.5 



0.07 
-0.31 
-0.55 
-0.90 



399 
399 
399 
399 



12.1 ±0.6 
10.6 ±0.6 
10.8 ±0.6 
11.5 ±0.6 



non- 


1 


6200 


2567 


46.26 


2225 


9.47 


2167 


0.03 


2889 


BALQs 


2 


4400 


3168 


46.04 


2469 


9.17 


2876 


-0.31 


3350 




3 


3400 


3693 


45.89 


2950 


8.92 


3861 


-0.55 


3312 




4 


2200 


3187 


45.60 


4981 


8.51 


3732 


-0.91 


3084 



'^ The property median value is followed by the number of objects in the bin and the observed fraction of BALQs from the total quasar 
population (in per cent). 
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§ 0.7 




fshift 



Figure 9. A zoom in on the Civ BAL profile for the compos- 
ite ratios presented in Figs. [6] and [8] The absorption profile be- 
comes broader and shifts to higher velocities as the Hen EW 
becomes lower, and the absorption becomes deeper at tishift > 
— 10,000 km s~^ as ayvi becomes redder. 



mainly in objects with a strong and b road C iv absorption 
(i.e., a high BI; see also Appendix K )) . iTrump et al.l (|2006l ') 
report that BALQs have broader emission lines than non- 
BALQs, which is consistent with the high Mgll FWHM bin 
having the largest /balQs (Table [6]). The small range in lu- 
minosity of the samples analysed here, prev ents detection 
of the Baldwin effect fe.g.. lGibson et al.ll2009l ). In this small 
range, we find a monotonic growth of /balQs with I/(3000A) 
(from 7.4 ± 0.4 to 15 . 2 ± 0. 8 per cent), which is also re- 
ported bv lAllen et al.1 ll201ll') and attrib uted to a S/N se- 
lection effect (see also lKnigge et al.ll2008 1. The trend found 
here is also affected by this effect, since there is a tight cor- 
relation in the low-z sample between the spectral S/N and 
I/(3000A), with the Spearman rank-order correlation coeffi- 
cient of rs ~ 0.57 and 0.66 for the BALQ and non-BALQ 
samples, respectively (null probability of logPr < —130). 
The median S/N per resolution element after smoothing 
(Section[2]) monotonically increases from < 40 for the lowest- 
L bin objects to > 100 for the highest-L bin objects, and 
some of the absorption is shallow enough to be affected by 
the S/N level (see Section [43)) . 



5 A PHYSICAL INTERPRETATION 

We describe below a possible physical context for our find- 
ings. First, we suggest how the Hen EW, which may be an 
indicator of the relative strength of the ionizing continuum 
beyond 4 Ryd, controls the Civ BAL iighift and /balQs (Sec- 
tion [S7T]). Then, we put forward a possible explanation on 
how auvi, which may be a reddening indicator, controls the 
C IV BAL depth and /balQs fSection [5^ . Finally, we briefly 
discuss connections to earlier studies (Section 15. 3|) . 



5.1 C IV BAL trends with the Hell EW 

What produces the strong trend between the Civ BAL iighift 
and the Hen EW? The Hen A1640 emission line is a re- 
combination line, and provides a measure of the number of 
Hell ionizing photons, i.e. photons above 54 eV, which are 
absorbed by the BLR gas. The Hell EW serves as a mea- 
sure of the relative strength of the EUV continuum above 
54 eV, compared to the near UV continuum. As the BAL 
outflow increases its radial velocity, its density must drop, 
and its ionization parameter must rise, potentially leading 
to complete ionization. Complete ionization halts further ac- 
celeration of a radiation pressure driven outflow. A higher 
Hell EW implies a higher ionization parameter at a given 
outflow density, and thus overionization of the outflow al- 
ready at a higher density, i.e. at an earlier stage of the wind 
acceleration, where the outflow velocity is smaller. Even if 
the outflow is not radiation pressure driven (say it is mag- 
netically driven), a harder ionizing continuum will overion- 
ize the outflow already at lower v elocities, but due to dif- 
ferent qualitative reasons (|Fukumura et al. 2010a b ) . In th e 
magnetically driven wind model of lFukumura et al.l (|2010al ^ , 
the wind density and velocity go as 1/r and l/\/r, respec- 
tively, where r is the distance from the ionizing source. The 
density law implies that the ionization parameter also goes 
as 1/r. Thus, in contrast with a radiation driven outflow, 
here the ionization parameter increases with increasing den- 
sity and decreasing r i.e., increasing velocity. However, for a 
harder ionizing continuum, the magnetically driven outflow 
gets overionized up to a larger r, i.e. down to a lower veloc- 
ity. Thus, the C iv BAL extends to lower maximal velocities, 
as expected for a radiation driven outflow, and as found in 
this study. 

To prevent overionization, iMurrav et al.l ([1995|) sug- 
gested the presence of a shielding gas, which blocks highly 
ionizing photons directed towards the outflow, preventing 
its overionization. Thus, ions such as C iv remain, despite 
the rise in the ionization parameter with increasing velocity. 
The observed relationship between Ushift and the Hell EW 
suggests that an outflow overionization does take place when 
the continuum which illuminates the BLR is hard enough. 
A large Wshift is obtained when the AGN ionizing continuum 
is soft, possibly enough to prevent overionization. Thus, one 
may not need to invoke a separate shield to prevent overi- 
onization. The suggestion that the observed weakness of the 
He II EW is sufficient to prevent overionization of the out- 
flow, can be explored quantitatively through photoionization 
modelling of wind model. 

The rise in /balQs with decreasing He ll EW is expected 
as a larger outflow volume is not overionized, leading to a 
larger CF of the BAL absorber. A magnetohydrodynamic 
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«shift [10'* km s^i 
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Figure 10. The dependence of the C IV BAL profile on ouvi s-iid the He II EW. The BALQ and non-BALQ samples are first binned into 
auvi bins (same as in Fig. [Sjl, and then each ouvi bin is binned based on the Hell EW. The upper left corner in each panel provides 
the median ayvi values for the BALQ/non-BALQ bins used to generate the plotted ratio. The corresponding values for the four Hell 
EW bins (in A) are indicated in the lower left corner of each panel. Note that all panels have the same y-scale. The C IV BAL profiles 
become shallower as ayvi becomes bluer, and at each ayvi bin tighift increases as the Hell EW decreases. The span in tighift increases 
as auvi becomes bluer. 



wind model may also produce a similar rise in the CF, as 
a softer ionizing continuum implies absorbing wind which 
extends to a smaller r, an d may produce a larger CF of 
the continuum source (JFukumura et al.ll2010bl ). It should be 
noted that since we are analysing average absorption pro- 
files, the measured CF has an inherent degeneracy between 
the CF as seen by an observer, the global covering factor of 
BAL gas (i.e., the CF as seen by the continuum source) and 
the fraction of B ALQs absorbing at a given iighift • 

Alternatively, the Hen EW may indicate the fraction 
of the sky of the ionizing continuum covered by the BLR 
gas fl instead of a measure of the strength of the EUV con- 
tinuum above 54 eV. However, there are various relations 
favouring the interpretation of the He ll EW a s an ionizing 
SEP i ndicator. The EUV slope is similar to Oox (|Zhene: et al.l 
I1997I : iLaor et all Il997l : iTelfer et al.ll2002l). The Opx is cor- 
related with luminosit y ( Strateva et al.l l2005l : ISteffen et al. 
I2OO6I : I Just et al.llioOTi '). and luminosity only jStern fc Laor 



l2012h . Simil arly, the He 11 A4686 EW is correlated with lumi- 
nosity only ([Boroson fc Greenll 19921 ). Thus, the ionizing con- 
tinuum gets softer with increasing luminosity, leading to the 
drop in the Hen EW with increasing luminosity. Also, lower 
lumin osit y objects have a flatt er EUV slopes ([Scott et al.l 
I2OO4I : cf. iTelfer et al.ll2002l and IShull et al.ll2012l ). In addT 
tion, a BLR il trend with luminosity is not expected to di- 
rectly affect Wshift- Future observations can explore directly 
whether the Hell EW is correlated with the EUV slope. 

In contrast with the Civ BAL, the Nv and O vi BALs 
do not show a prominent difference in their absorption pro- 
files with the Hen EW (Fig. (5]). The blue wing of Ovi may 
be affected by Nlll A991 absorption, and to a lesser extent 
by C III A997 absorption, while the N v absorption is affected 
by Ly a absorption. Note also that the N v and O vi ions are 
produced and destroyed by photons at energies significantly 
above 54 eV, in contrast with the C iv ion, where these en- 
ergies (47.9-64.5 eV, Section [3| are just below and above 



© 2013 RAS, MNRAS 000.[TH22l 



16 A. Baskin, A. Laor and F. Hamann 



g-<0.9 



g 0.1 



ttj 




(5 



0.7 



4"> (low) 
Mg 11 FWHM bins 



A/bh bins 

-2 -1 

''shift [10* kin s~^ 



i{3000A) bii 




1^' (liigh) 



-2 -1 

J'shitt [10" km s~i] 



Figure 11. Same as Fig. |9] but the samples are binned based 
on the Mgll FWHM, L(3000A), A/bh and L/LEdd- Since the 
L(3000A) range is small (see TableO, the trend in the Civ ab- 
sorption with Mbh and L/LEdd is a simple outcome of the corre- 
lation between the C IV absorption and the Mg II FWHM [panel 
(a)]. Note that all presented trends are weaker than the trend 
between the Civ BAL and the Hell EW and ctyvi (Fig-[21l- 



4 Ryd. One thus needs to further explore with photoion- 
ization calculations the relation between the Hell emission 
EW, and the expected N v and O vi columns relative to the 
C IV column, as a function of the ionizing spectral shape. 



5.2 C IV BAL trends with auvi 

What produces the trend between the absorption depth and 
Quvi ? The value of auvi may be interpreted as a viewing 
angle indicator, if it is affected by reddening, and if the dust 
tends to reside in the symmetry plane of the system. A sup- 
port for this scenario is provided by a relation between the 
optical-UV slope (qp.uv) and the degre e of white light polar- 
izatio n, found by Baskin & Laor ( 2005) in the complete sam- 
ple of lBoroson fc Green (1992 ) PG quasars. A redder Qo,uv 
is associated with a higher polariza tion, as expecte d for a 
system observed closer to edge on. iBaskin fc Laorl (|2005l ) 
also show that the change in Qq uv is consistent with dust 
reddening (see also lStern fc Laorll2012l . fig. 17 there). In this 
interpretation, if there is a clumpy and planar distribution of 
absorbers, then a redder auv corresponds to viewing angles 
closer to edge-on, for which there is a larger probability that 
our line of sight intersects a UV absorbing gas cloud, and 
also a dusty gas cloud, i.e. a larger BAL CF and a larger 
reddening. This explains both the deeper absorption at a 
given Ushitt, as a larger fraction of the central source is ob- 
scured, and the rise in /balQs, as a larger fraction of lines 
of sight passes through an absorber. It is interesting to note 
that the median absorption depth increases from ~ 5 per 
cent at the shallowest auvi, to ~ 20 per cent at the steepest 
Quvi (Fig.O, while /balQs changes similarly from 6.9 ±0.4 
to 21 ± 1 per cent with auvi (Tabled. 

BALQs observed at a larger inclination are also likely 
to have a larger absorbing column. However, since all lines 
are most likely saturated, they are not expected to show 



a rising absorbing EW. However, X-ray observations may 
reveal larger absorbing columns in steeper auvi BALQs. 

The inclination interpretation of auvi implies that 
BALQs and non-BALQs, at a given auvi, are observed at 
the same inclination. What separates BALQs from non- 
BALQs is then whether an outflow component happens to lie 
along our line of sight. This scenario requires an azimuthally 
asymmetric outflow structure, whic h is also hinted at by 
the a bsorption profile variability (e.g.. lCapellupo et al.ll201ll . 
|2012| ). 

The inclination interpretation of auvi can also explain 
the lower Wshift observed in redder BALQs (Figure fTO|l . View- 
ing a BALQ close to edge-on increases the probability that 
the line of sight goes through the base of outflowing ab- 
sorbing gas, where the projected velocities are lower. A 
trend indeed obs e rved in recent disc wind simulation by 
iGiustini fc Progal (|2012l . figure 5 there). In contrast, the 
«shift > 10,000 km s~* part of the outflow shows only a 
weak auvi dependence (Figure |9]), which suggests the higher 
velocity outflow component is more spatially extended, and 
thus less inclination dependent. 

We note in passing that binning by the shorter wave- 
lengths auvs leads to C iv profiles which are very similar to 
those derived from the Hell EW binning, which is in con- 
trast with the result for the auvi binning. Thus, auvs is 
not orthogonal to the Hell EW, as auvi is. This can also 
be seen in Table |3l where the two auvs bins show a fac- 
tor of 4 difference in their median Hell EW, in contrast 
with their similar values in the auvi bins (Table |4}. Thus, 
auvs appears to be related to the ionizing SED, as measured 
by the He ll EW, while auvi is independent of the ionizing 
SED, and likely provides a measure of the dust extinction 
(and indirectly the inclination). This qualitative difference 
is not unexpected, if the optical-UV continuum is produced 
by accretion disc emission. At long wavelengths, accretions 
discs are expected to show a universal SED slope, as we are 



observing the universal slope L^ oc 



,1/3 



part of the disc 



emission (realistic model slopes are not that flat), which is 
independent of Mbh and L/LEdd- At short enough wave- 
lengths, the slope starts to probe the position of the accre- 
tion disc spectral turn over re gion, and is thus a measure of 
the ionizing SED (e.g. iDavis fc Laor., 2011. ') . 



5.3 The luminosity dependence of Vmax 

iBrandt et all (|2000l) and iLaor fc Braiidtl (|2002| ) found a 
clear correlation of the maximal outflow velocity Umax and 
the AGN luminosity in soft X-ray weak objects. This re- 
lation is quite robust as no Seyfert-level AGN is found to 
reach Umax ~ 10,000 km s~*, which is common in BALQs. 
A pure luminosity tre nd is simply explain ed by radiation 
pressure driven winds (jLaor fc Brandtll2002l . sec. 3.4 there), 
which gives Umax oc L^' * for a wind launched from the BLR 
with a constant force multiplier. The observed steeper rela- 
tion Umax oc 7^062±o.08 ^^^ |^g interpr eted as an indication o f 
a force multiplier that rises with L (jLaor fc Brandtll2002l ). 
Such a rise may be caused by a softer ionizing SED at a 
higher L, which produces less overionization, and thus a rise 
in the force multiplier. A softer ionizing SED at a higher L 
is consistent with the observed inverse relation between the 
Hell EW and L. The range in AGN luminosity is too small 
here to test the luminosity dependence of Umax (but see ten- 
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tative evidence in Appendix [Bj . However, the trend found 
here for Vshift with the Hen EW, and the strong luminosity 
dependence of the He ll EW, may indicate that the observed 
relation between Dmax and the AGN luminosity found by 
iLaor fc Brandj ([2002^, may be partly driven by the Vshitt 
versus Hen E W relationship found h ere. A quick inspection 
of the data in lLaor fc Brandtl l|2002l ) reveals that at a given 
luminosity, the highest Wmax is reached by the lowest Hen 
EW objects (which are also the soft X-ray weak quasars). 
However, a large sample which covers a large range in lumi- 
nosity and the Hen EW is required to clearly separate out 
the dependence of Ushift on L and on the He ll EW. 

Inspection of Fig. [6] reveals that the Civ emission 
profiles get weaker and blueshifte d with decreasing He ll 
EW, consistent with the finding of iReichard et al.l ([2003) • 
[Richards et al.l (|201ll . fig. fl there) demonstrate the tight 
relation between the He ll EW and the C iv emission 
EW and asyrametry, which they interpret as an indica- 
tion for_einissioii_^x)m a wind component in the BLR (see 
also lLeighlv fc Moorell2004l : iKruczek et al.ll201ll ). Above we 
found that the BAL wind component becomes more promi- 
nent as the EUV ionizing continuum gets weaker, which sup- 
presses the overionization of the wind. The suggested BLR 
wind component may very well be related to the base of 
the BAL outflow, which is most likely fed by the BLR gas. 
As the wind accelerates from the BLR, its density and thus 
emissivity drops, but it remains visible as a BAL through 
its resonance line absorption. A weaker Hen EW implies 
less overionization, which allows a higher wind emissivity 
closer to the base of the wind, producing a blushifted C iv 
emission component. Unlike the blue wing of the C iv emis- 
sion that is probably produced by a matter-bounded wind 
component, the red wing of CiV is probably prod uced by 
a non-wind (disc) component (JRichards et al.l 120111 ). which 
is likely radiation bounded. As the relative strength of the 
EUV ionizing continuum becomes weaker, there is less pro- 
duction of C^^ in the radiation-bounded disc component, 
which leads to a weaker emission in the C iv red wing with 
decreasing Hen EW, as observed. 



6 CAN DUST EXTINCTION EXPLAIN THE 
REDDENING? 

The median BALQ SED is redder than the non-BALQ SED 
(Fig. [l|. The difference between the two SEDs appears to 
be caused by the overall spectral slope, rather than local 
features. Is the difference in spectral slopes intrinsic to the 
illuminating source, or is it a result of extinction by fore- 
ground dust which is more common in BALQs? 

Figure [12] presents corrections of the BALQ composite 
by several possible dust extinction laws. The BALQ compos- 
ite is corrected to match the non-BALQ composite by adopt- 
ing three types of dust extinctio n laws: Milky Way (MW; 
ICardelli. Clavton fc Mathid [19891 ') . Large Magellanic Cloud 
(LMC) and Small Magellanic Cloud (SMC; the latter two 
laws are from .Gordon et al.ii2003i ')l''l The upper panel shows 



° All three extinction laws are reported only down to Arest = 
1000 A. Here we assume the analytic functions hold down to 
Arcst = 800 A. Given the low S/N of the R spectrum at Arcst < 
1000 A, the exact form of the extinction law used is not important. 



the de-reddening for the high-z sample. One cannot discrim- 
inate between the three extinction curves. The lower panel 
shows the de-reddening for the low-z sample. Here one can 
clearly exclude a MW- and LMC-like extinction laws by ~ 13 
and 9 standard errors in the median R at Arest ~ 2200 A, 
as they produce a strong bump at 2000 < Arest < 2400 A, 
which is not observed [j 

A preference for SMC dust was found in various earlier 
studies of AGN. ISpravberrv fc Folta l|l992t ) reach a simi- 
lar conclus ion for th e reddening of LoBALQs compared to 
HiBALQs. IChkman et al. (2012 ) find that the extinction of 
dust-reddened quasars in the Arest= 4000 A- 2.4 ^m range is 
best described by an SM C-like law (see also a brief overview 
of additional studies by IPitman. Clavton fc Gordon! [2000|) . 
iHopkins et al.l (^2004) report that the reddening of quasars in 
general is consistent with an SMC-lik e extinction law , while 
the MW and LMC laws are excluded. [Richards et al] (120031 ) 
also successfully fit dust-reddened non-BALQs assuming an 
SMC extinction law. Thus, the result derived here is clearly 
not new. However, these earlier studies are all based on the 
observed SED. The observed SED happens to show a shal- 
low and broad dip at A ~ 2200 A, produce by adjacent line 
emission, which can mimic some effect of the 2200 A ex- 
tinction bump. The ratio plots used here eliminate all the 
intrinsic spectral features, and thus provide a strong con- 
straint (to a level of ~ 1 per cent) on the possible extinction 
curve shape associated with the BAL outflow. 

The SMC best-fit Av is evaluated independently for 
the high- and low-z samples. Both samples yield Ay ~ 
0.06 ± 0.01 mag, where the range implies a deviation that 
is consistent with the standard error in the median R spec- 
trum. The best fit Ay is estimated by eye inspection, re- 
quiring the non-BALQ and dereddened BALQ composites 
to match in wavelength regions unaffected by strong line 
emission and absorption, i.e. Arcst ~ 1080 and 1800 A for the 
high- and low- 2; samples, respectively. The Ay = 0.06 mag 
reddening fits the top three Hen EW bins of the low-2 
sample (Fig. O, while the lowest Hen EW bin requires 
Av = 0.09 mag. 

The Ay values found here are mostly consistent with 
what was previously reported for BALQs (^Reichard et al 



[2003; Maddox ct al. 2008; Gibson et al. 2009). Alle n et al 
(|201H ) find a trend between Ay and z. The Ay derived 
here i s consistent with the value evaluated by JAllen et al.l 
(J201H ) for 1.6 < z < 2.0 range (our low-z sample), but it is 
~ 3 times smaller than the value for the 3.8 < z < 4.5 range 
(our high-z: sample). We suspect the 3.8 < z < 4.5 result of 
JAUen et al.l (|201ll ) is an overestimate, as it over-corrects the 
BALQ Nv line emission (see fig. 22 there). 

The implied H column is A''(H) « 10^^_cm~Vusing the 
iV(H) = 5.8 X 10'^^ Ay /Ry relation (Drainc"201l'), where 
Ry = 2.74 for the SMC extinction law (Gordon et al. 2003i ). 
Given the above constraint on the H° column, based on the 
Lyman edge for a CF — 1 absorber (Section 13. ip . the H 



^ The standard error in the median R is evaluated as follows. 
First, the standard deviation of fx at 1700—1720 A is calculated 
for the median BALQ and non-BALQ spectra. Then, the error as 
a function of Arest is estimated by assuming the measured error 
is dominated by photon statistics i.e., the standard deviation is 
scaled by a//a (Arest)- Finally, the error in the R spectrum is 
calculated by the standard error progression. 
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Figure 12. High- and low-z BALQ composites corrected for dust extinction assuming different extinction laws. The \ow- z comp osite 
is normaUzed at 3000 A for presentation purp oses. Three ex t inctio n-laws are adopted: MW {Rv = 3.2) from lCardelli et al.l 1 119891 ). and 
LMC {Rv = 3.4) and SMC {Rv = 2.74) from lCordon et al.l | |2003| '|. The best fit Ay is estimated by eye inspection (see text). The error 
bars indicate the standard error in the median value. For the high-2 sample (upper panel), all extinction laws are possible, and the MW 
and LMC extinction corrections are essentially identical. However, the low-z: R spectrum (bottom panel) excludes the MW and LMC 
extinction laws at a high significance level, as they produce a strong bump at 2000 < Arest < 2400 A, which is not observed. Note that 
both the low- and the high-2 composites are consistent with an SMC-likc dust extinction with the same Av value of 0.06 mag. 



column associated with the dust must be nearly fully ion- 
ized. If the gas is diffuse, the ionization can be maintained 
by the AGN on the host- galaxy scale. Since non-BALQs are 
also d ust reddened (e.g.. [Richards et al.ll2003l : IStern fc Laon 
[2013), the above estimate gives only the excess extinction in 
BALQs. To get the absolute extinction value, we evaluate 
the low-z BALQ composite reddening relative to the bluest 
non-BALQ quartile composite, which is more likely to rep- 
resent the intrinsic non-reddened SED. This yields Av = 
0.12 mag and A'^(H) ~ 3 x 10^" cm~^. A similar analysis for 
the high-z sample, utilizing the bluest non-BALQ compos- 
ite, yields ^v = 0.10 mag and iV(H) Ri 2 x 10^" cm"^ 

Since the SMC extinction curve is featureless, an intrin- 
sically redder SED for BALQs cannot be excluded. The only 



support for the reddening interpretation is the remarkably 
good match, i.e., a ratio of unity for the dereddened spectra 
throughout the observed range, from 1000 to 3000 A. 



7 CONCLUSIONS 

We analyse the average absorption properties of high- and 
low-z samples of BALQs from the SDSS DR7, which cover 
the wavelength range of 800-3000 A. The absorption proper- 
ties are derived by taking the ratios of BALQ median spectra 
to the median spectra of matched samples of non-BALQs. 
We explore the absorption properties for a number of sub- 



© 2013 RAS, MNRAS 000.[TH22l 



Average BALQ absorption properties 19 



samples selected based on various emission properties. We 
find the following. 

(i) No Lyman edge associated with the BAL absorbing 
gas is detected (r < 0.1). Thus, on average, CF < 0.1 for a 
partially ionized absorber in BALQs. 

(ii) The average absorption EW increases with the ion- 
ization state, from an absorption EW < lA for C ii A1335, 
to 5.8 for Si IV, 9.3 for Civ, 24.6 for Nv and 25.2 A for O vi. 
This may indicate a rise in the covering factor of the BAL 
outflow with a rise in the ionization state. 

(iii) The Hell emission EW controls the typical «shift of 
Civ BAL, which increases from v^mtt < 7000 km s^^ for 
EW > 6 A to Dshift > 15000 km s~^ for EW < 1 A. The 
Hell EW does not affect the absorption depth. The Hell 
EW may indicate the strength of the EUV. A lower Hell 
EW then implies a lower ionization of the outflow, which 
allows the outflow to reach higher velocities before being 
overionized. One may not need to invoke a shield to pre- 
vent over-ionization of the outflow, as over-ionization may 
be taking place in high He ll EW objects, while in low He ll 
EW objects a shield may not be needed. 

(iv) The value of auvi controls the C iv peak absorption 
depth, which increases from ~ 0.15 for Quvi ~ —0.5 to 
~ 0.45 for auvi « -1.5, at WsMft < 10,000 km s~\ The 
value of auvi does not affect the typical iishift and the ab- 
sorption profiles at Ushift > 10,000 km s~^. The value of 
auvi niay control the average inclination angle of the sys- 
tem, and thus the covering factor and absorbing column for 
a planar outfiow, in particular at the lower velocities, which 
likely originate close to the base of the outflow. 

(v) The Hell EW and auvi also control the fraction of 
AGN which are BALQs. The fraction rises from 4.2 ± 0.4 
per cent in blue (auvi ~ —0.5) and strong Hell (EW « 
6 A) AGN, to 31 ± 4 per cent in red (auvi ~ —1.5) and 
weak He II (EW < 1 A) AGN. Also, most of the LoBALQs 
(~ 75 per cent) have Hell EW < 1 A. The Hell EW may 
control the global CF of the BAL outfiow, as seen by the 
ionizing source, as a lower Hell EW allows the outflows to 
extend to larger scales before being overionized. On the other 
hand, auvi niay control the CF of absorbers along the line 
of sight, which increases with increasing inclination and the 
associated reddening. 

(vi) The median SED of BALQs is consistent with ex- 
cess reddening of purely SMC dust with Av = 0.06 mag 
compared to non-BALQs, or a reddening oi Av =0.12 mag 
compared to the bluest non-BALQs. Given the lack of as- 
sociated Hi absorption, the dust is embedded in an ionized 
medium, possibly on the host-galaxy scale. LMC and MW 
dust are excluded by ~ 9 and 13 standard deviation of the 
median value, respectively. 

The above interpretation of auvi as a viewing angle in- 
dicator can be tested by looking for a trend between auvi 
and the continuum polarization, in both BALQs and non- 
BALQs. One can also look for the expected relationship be- 
tween auvi and the X-ray absorbing column in BALQs. The 
relationships of VsMit and /balQs with the Hell EW, to- 
gether with the inverse relationship of the Hell EW with 
L, imply a steeper than linear rise in the kinetic wind lu- 
minosity with L, which may be relevant for feedback in the 
highest L AGN. 
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APPENDIX A: THE LoBALQ R SPECTRUM 

As a sidenote, we use this opportunity to briefly present the 
LoBALQ median absorption properties based on the low- 
2 BALQ sample. The LoBALQ R spectrum is derived by 
forming a matched non-BALQ sample, which is matched by 
the Hen EW distribution. Although most of the LoBALQs 
have small values of Hell EW (the whole sample median 
EW is —1.8 A), several LoBALQs (13 out of 56) have some- 
what larger Hell EW (> 1 A; Tabled, and the matching 
method used for HiBALQs (i.e., matching in the Hell EW 
range) cannot be utilized here. To produce LoBALQ and 
non-BALQ composites with a matched Hell EW distribu- 
tion, we carry out the following procedure. 

(i) The LoBALQ Hell EW distribution is calculated us- 
ing 2 A bins (EWbin). At each EWbin, the number of LoB- 
ALQs with He II EW within the EWbin range is divided by 
the total number of LoBALQs. 

(ii) The non-BALQs are assigned to EWbin based on their 
Hen EW value. This yields the measured absolute non- 
BALQ distribution, iVmeas (EWbin). 

(iii) The LoBALQ distribution is multiplied by the total 
number of non-BALQs. This yields the desired absolute non- 
BALQ distribution, iVdos (EWbin). 

(iv) The non-BALQs are drawn for each EWbin based 
on the following three criteria. If A'mcas (EWbin) = 
A'dcs (EWbin), then all non-BALQs in that EWbin are drawn. 
If A^moas (EWbin) < A^dos (EWbin) , then all non-BALQs are 
randomly drawn several times until the number of drawn ob- 
jects equals A"dos(EWbin). If A^moas (EWbin) > A"dos(EWbin), 

then only A'^doa (EWbin) of non-BALQs is randomly drawn. 

(v) The median spectrum of the drawn non-BALQs is 
evaluated. 

(vi) The resulting median non-BALQ spectrum depends 
on the particular non-BALQs that are randomly drawn. To 
eliminate this dependence, the procedure in items (iv) and 
(v) is repeated 300 times. A median of the 300 median spec- 
tra is adopted as the matched non-BALQ composite. 

The above procedure produces a matched non-BALQ sam- 
ple with a median Hell EW of —1.3 A (compared to —1.8 A 
for LoBALQs). The LoBALQ and the matched non-BALQ 
composites have auvi of —1.99 and —1.38, respectively. The 
resulting Civ BAL EW, v'^^f^ and a are 33.9 A, -10,500 
and 5300 km s^^, respectively. 

Figure [ATI presents the LoBALQ R spectrum corrected 
for reddening (see below), where we also indicate the labora- 
tory wavelength of Mgii AA2796.4, 2803.5, Aim AA1854.7, 
1862.8 and several Fell absorption lines following iRichardj 
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l|200ll . table 3 there). Note the presence of prominent Fell 
BALs in the LoBALQ R spectrum (absorption EW > 7 A) . 
The LoBALQ C iv BAL profile is consistent with the re- 
lationships, which are found for the low-z HiBALQ sam- 
ple, between f shift and the Hell EW, and between Civ 
peak absorption and auvi (Sections 14.11 and 14.21 respec- 
tively). The LoBALQ composite has a redder quvi and a 
deeper C iv peak absorption than the reddest HiBALQ bin 
(cKuvi = —1.99 versus —1.51; compare to Fig. [S]). The Civ 
BAL Ushift of LoBALQ composite is similar to Ushift of the 
lowest Hen EW HiBALQ bin (Fig. EH), ahhough the LoB- 
ALQ median Hell EW is lower (0.9 versus -1.8 A). The 
similarity in Ushift may be a result of the softer ionizing SED 
(i.e., lower Hell EW for LoBALQs) being 'too soft' to pro- 
duce a C IV absorber at higher iighift . The higher VsMit ab- 
sorber is probably found at larger distances from the ioniz- 
ing source than a lower Wshift absorber, and the soft ionizing 
SED at those distances results in an ionization parameter 
that is too low to produce C iv absorption. 

A comparison between LoBALQ and HiBALQ compos- 
ites supports the scenario of a smooth transition between 
HiBALQs and LoBALQs, where the latter mostly occupy 
the reddest auvi and the lowest Hell EW parameter-space 
(Section l4.3p . Figure lATI presents the bluest and highest Hell 
EW composite, the lowest He ll EW composite (Fig. [6| cor- 
rected for reddening (Section[6l), and the reddest and lowest 
Hen EW composite. The first composite has EW = 6.1 A 
and auvi = —0.48 (Table [5]), and its R spectrum has no 
detectable Aim and Mgll BALs. The dip at Mgll labora- 
tory wavelength results from a mismatch in Mg ll peak line 
emission between the matched BALQs and non-BALQs, and 
likely does not indicate true absorption. The second compos- 
ite has EW = 0.9 A and Quvi = -0.96 (Table [g]), and has 
a marginal Aim BAL and non-detectable Mgll BAL. The 
third composite has a Hell EW — —0.6 A and auvi = —1.51 
(Table [SJ, and its R spectrum has a prominent Aim BAL 
and a marginal Mgll BAL. There are additional evidence 
of Aim BAL strength dependence on both Hell EW and 
auvi- HiBALQ binning based on Hell EW (cvuvi) produces 
composites with similar auvi (Hell EW), but with different 
Aim BAL strength (Figs. [6] and [8|. Note that the R spec- 
trum in the Mgll region in Figs. [B] and [S] is dominated by a 
mismatch in the Mgll peak line-emission, which hinders a 
detection of any weak Mgll BAL. The Aim BAL appears 
first, as the ionization parameter becomes lower, since Mgll 
ions are destroyed by 15 eV photons, while Aim ions are pro- 
duced by 19 eV photons. As the HiBALQ properties (i.e.. 
He II EW and auvi) start to resemble those of LoBALQs, low 
ionization BALs begin to appear in the HiBALQ spectrum, 
and it becomes more similar to the LoBALQ spectrum. This 
suggests LoBALQs are the extreme high inclination and soft 
ionizing SED members of the BALQ population, rather than 
a distinct population of AGN. 

The reddening of the LoBALQ R spectrum is consis- 
tent with an SMC dust extinction oi Ay = 0.22 mag i.e., 
E(B- V) = 0.08. This reddening is sim ilar to E(B- V) = 
0.077 reported bv lReichard et al.l (|2003t ) for the SDSS early 
DR, but is s malle r than E(B—V) = 0.14 reported by 
ICibson et al.l ((20091) for the DR5, where in both studies the 
reddening is derived by comparing the LoBALQ composite 
to a composite of the whole non-BALQ sample. If the whole 
non-BALQ sample, rather than the Hell EW matched one. 



is utilized to calculate the LoBALQ R spectrum, then the 
ICibson et al.l ([20091) result is approximately recovered i.e., 
Av = 0.41 mag. 



APPENDIX B: TENTATIVE EVIDENCE FOR 
RADIATION PRESSURE DRIVING 

Figure IBll presents tentative evidence that the C iv BAL 
absorber is driven by radiation pressure. The low-z sam- 
ple is divided into 4x4 bins based on the Hell EW and 
L(3000A). The BALQ sample is first divided into four Hell 
EW 'parent' bins, and matched non-BALQ bins are con- 
structed (same procedure as in Fig.[B|. Then, each Hell EW 
parent bin is divided into four L(3000A) bins, with matching 
non-BALQs from the corresponding Hell EW non-BALQ 
bin. Each BALQ bin contains by construction equal number 
of objects (~ 100) for all 'parent' bins. If the BAL outflow 
is driven by radiation pressure, then Ushift is an increasing 
function of the force multiplier M and luminosity L, as- 
suming th e outflow launching r adius is mainly a function 
of L (e.g., iLaor fc Brandtl [2003 : see also Section [5|. M is 
mostly set by the SED, and keeping the Hell EW constant 
should also keep M approximately constant i.e., each panel 
of Fig. IBll groups objects with a similar M. The range in 
I/(3000A) of our sample is too small (< 1 dex) to produce 
highly significant differences in Vshift between L bins for a 
given Hell EW (i.e., M). However, note that for the weakest 
two Hell EW bins (Fig. IBll two lower panels), the highest 
L bin reaches larger Wshift than the other three bins. This 
larger VsMit for a larger L is expecte d, if the BAL outflow 
is dr iven by radiation pressure (e.g., [Proga. Stone fc Drewl 
Il998h . 

The highest Hell EW bin shows very similar iishift 
as a function of L. Here the peak absorption occurs at 
^ahift ~ 3000 km s~^, and the absorbing material may still 
be confined close to the BLR, i.e. to the base of the out- 
flow. Overionization occurs just as the outflow starts to build 
up, before radiation pressure has time to build up a signif- 
icant outflow component, beyond the initial dispersion in 
the velocity. Only when the Hell EW is low enough, to al- 
low signiflcant acceleration, the effect of L starts to appear. 
Clearly, large sample s which span a large range in L (e.g. 
[Laor fc Brandtll2002l ). are needed to explore the radiation 
pressure interpretation, and clearly separate out the effect 
of L and the He ll EW on the outflow. 
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Figure Bl. The dependence of the Civ BAL profile on the Hell 
EW and L{3000A). The BALQ and non-BALQ samples are first 
binned into Hell EW bins (same as in Fig. |6)l, and then each 
Hell EW bin is binned based on L(3000A). The upper left cor- 
ner in each panel provides the median Hell EW (in A) for the 
BALQ/non-BALQ bins used to generate the plotted ratio. The 
corresponding values for the four IogL(3000A) bins (in erg s~^) 
are indicated in the lower left corner of each panel. Note that 
for the lower two Hell EW bins (two lower panels), the highest 
L(3000A) composite (blue line) has the deepest absorption trough 
at Dshift ^ -2 X 10"' km s'^. 
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